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Abstract 



Within the nuclear Fermi-liquid drop model, quantum and thermal fluc- 
tuations are considered by use of the Landau- Vlasov-Langevin equation. The 
spectral correlation function of the nuclear surface fluctuations is evaluated 
in a simple model of an incompressible and irrotational Fermi liquid. The de- 
pendence of the spectral correlation function on the dynamical Fermi-surface 
distortion is established. The temperature at which the eigenvibrations be- 
come overdamped is calculated. It is shown that, for realistic values of the 
relaxation time parameter and in the high temperature regime, there is a 
particular eigenmode of the Fermi liquid drop where the restoring force is 
exclusively due to the dynamical Fermi-surface distortion. 
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I. INTRODUCTION 



The collective dynamics and the dissipative properties of the nuclear Fermi liquid depend 
in many aspects on the dynamical distortion of the Fermi surface in momentum space. As 
is well-known, taking into account this distortion allows the description of a new class of 
phenomena, the most famous of which are the giant multipole resonances. Furthermore, the 
scattering of particles on the distorted Fermi surface leads to relaxation of collective motion 
and gives rise to nuclear viscosity [|I|. This so-called collisional mechanism of relaxation is 
influenced very strongly by the Pauli blocking in the scattering of nucleons close to the Fermi 
surface and depends on temperature and retardation effects PQj. We will consider below 
both of them, concentrating on the fact that the retardation effects are especially important 
for a proper description of the transition from the zero-sound regime to the first-sound 
regime in the excited (heated) nucleus. 

In the present paper, we are interested in the spectrum of fluctuations in shape variables. 
The precise form of such spectra can be expected to depend on the parameters of the model, 
such as the collision time, and, especially, on the memory effects. Here, we want to study 
these dependencies as one step to our ultimate goal of determining the model parameters 
from a comparison with experimental data, as might be possible due to a relation of the 
afore mentioned spectra to 7— spectra. 

In what follows, we combine the thermal and quantum fluctuations by means of the 
fluctuation-dissipation theorem. Such an approach presents a convenient connection be- 
tween different regimes of collective motion such as the quantum zero-sound regime at zero 
temperature and the collisional first-sound regime in a hot system. 

This paper is organized as follows. In section II we suggest a proof of the Langevin 
equation for nuclear shape variables, starting from the Landau- Vlasov kinetic equation. In 
the derivation, the main features of the dynamical distortion of the Fermi surface are taken 
into account. In section III, we discuss the collision integral, for which a form is adopted that 
includes memory effects. The limiting cases of the zero-sound and the first-sound regimes 
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are discussed in section IV. Results of numerical calculations are presented in section V. We 
conclude and summarize in section VI. The Appendix provides a derivation of the relevant 
correlation functions for the fluctuating forces in the case of non-Markovian processes based 
on the (second) fluctuation-dissipation theorem. 

II. SPECTRAL CORRELATION FUNCTION FOR SHAPE FLUCTUATIONS IN 

A FERMI-LIQUID DROP 

To consider fluctuations of the collective variables in a Fermi-liquid drop we start from the 
following kinetic equation for the small deviation 5n of the distribution function n = n(r, p; t) 
from the one in equilibrium, Ueq, taking into account a random force 

d 

—Sn + L6n = I[6n]+y. (2.1) 

On the left hand side of ( |2.1|) , the operator L represents the drift terms including the 
selfconsistent mean field U, 

m 

On the right hand side, I[6n] is the linearized collision integral and y = y{r, p; t) is a random 
variable representing the random force. As such, its ensemble average vanishes, (y) = 0, 
while its second moment can be related to properties of the collision term, as shown in the 
Appendix. 

To derive the equation of motion for the shape variables, we will follow the nuclear fluid 
dynamic approach of Ref. [Q, and take into account the dynamic Fermi surface distortion 
up to multipolarity 1 = 2: 

l,mi 

Here e is the quasiparticle energy |T|]. A generalization of our approach to the case of an 
arbitrary multipolarity / of the Fermi surface distortion can be done in a straightforward 



way, see Ref. 0. Using Eqs. (|2.1| ) and (|2.2|) we can derive a closed set of equations for 



the following moments of the distribution function, namely, local particle density p, velocity 
field and pressure tensor P^^, in the form 

d 

g^5p= -WuipegUu), (2.3) 
^Pe^Q-U, + Pe,'^u{j^)jp + V^P^^ = 0, (2.4) 

d 2 

— P^^ + Peq{VuU^ + Vf,U,y - -6uf,VaUa) = hfi + V^fi, (2.5) 

exploiting the conservation of particle number and momentum in collisions between the 
particles. Here S is the internal energy density, which is the sum of the kinetic energy 
density of the Fermi motion and the potential energy density associated with the nucleon- 
nucleon interaction. The equilibrium pressure of a Fermi gas. Peg, is given by 

^- = i/(2^^'"- ^'-'^ 

P^^ is the deviation of the pressure tensor from its isotropic part due to the Fermi surface 
distortion 

^!^f^ = ~~J ^^^|p(P^ -"^^'^)(Pa' -"^«M)(^)ggII'^^2m(r,t)F2„(p), (2.7) 
/j,^ is the second moment of the collision integral 

and Ui^fj^ gives the contribution from the random force 

1 f dp . . 

Using the Fourier transformation for the pressure 

K,it) = J ^e-*P;,,^ (2.10) 
and similarly for the other time dependent variables we find the solution to Eq. ( p.5|) as 



= '"^i^^JJy PeqK^.,u + Y^/j^^^ + ^^r)y,^^^, (2.11) 
where we used the symbol 
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Kii,u; = '^uXti,u + '^^lXv,u - ^5^p^\X\,u; (2.12) 

for this combination of gradients of the Fourier transform Xu,u) of the displacement field. The 
time derivative of x(r, t) is defined as the velocity field, hence 

Uy^uj = -iuJXu,uj- (2-13) 

To obtain Eq. ( |2.11| ) we have also used the fact that the tensor Jj,^, Eq. ( |2.8| ), can be 
reduced to 

T ^ 

due to our restriction to quadrupole deformation of the Fermi surface. This is because the 



/ = and 1 components of the expansion (2^) do not contribute to the collision integral, 



reflecting the conservation of particle number and momentum in a collision. Note that the 



form ( 2.14 ) is also correct for a Non-Markovian collision term, with the collision time r being 
dependent on the frequency uj. (For convenience we will omit this frequency dependence of 
r in our notations.) 

From Eqs. (|2.3| ), (|2.4|) and (p.ll|) we find the equation of motion for the displacement 
field ill the form 

- Peq(^'^Xu,uj + ^Xu,u; = V^((Tj,^,^^ + Sj,^,^), (2.15) 

where the conservative terms are abbreviated by 

^Xu,cv = -—peq'^u(Y^) V ^peqXf^,u - 'i^T^ ( ) V^ — A,.^,^;, (2.16) 

m ^dp^^eq •J \1 — iujtJ m 

is the viscosity tensor 

(yuf,,u = -i{uj/m)r]{uj)Ay^^^ (2.17) 



with the viscosity coefficient 



v{iu) = Re[- — Pe,, (2.18) 

,1 — tUOTj 



and Sufj,,uj is the random pressure tensor 

r(l + iur) 

The correlation properties of s,y^,(^ can be obtained for the general case where we also take 
into account retardation and memory effects in the system, see the Appendix for details. 
Using the correlation properties of the random tensor y^^^a;, we find for the ensemble average 
of 

the result, see Eq. (A. 17), 



Att 2 
—E{uj,T)r]{uj)6{r - r')6{uj + uj')[6^^>6f,f,' + 6^^:6f,^> - -6^f,6^:f,>], (2.20) 



where 



B(u,,r) = ^coth^. (2.21) 

We have pre j served the constant h in Eq. ( p.21| ) in order to stress the fact that both quantum 
and thermal fiuctuations are involved in Eq. ( |2.2CI| ) [^|[. 

To apply ( p.l5| ) to a finite system, we shall assume a sharp boundary of the Fermi liquid 



p = PoQiR{t)-r). (2.22) 



Below, we will also assume this simple form of the particle density p, Eq. (|2.22| ), at non- 
zero nuclear temperatures T 7^ 0. In a hot nucleus, the particle density parameter po in Eq. 
( p.22| ) is temperature-dependent, po = Po(T). However, the form (|2.22| ) ignores the existence 
of a nucleon vapour. 
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For the description of small amplitude oscillations of a certain multipolarity L of a liquid 
drop we specify the liquid surface as 

r = R{t) = Ro[l + E«iM(t)>LM(^,0)]. (2.23) 

M 

We write the displacement field X;^(r, t) for an incompressible and irrotational flow, V^Xu = 
0, as i 

Xi.(r, t) = J2 aLMA^)^LMit), (2.24) 

M 

where 

(^lmA^) = 74:^2 V.(r^riM(^, </>)). (2.25) 
LKq 

Multiplying Eq. (|]l|) by ma^^j^, summing over u and integrating over r-space, we obtain 
the Langevin equation for the collective variables, 

- Uj'^BLaLM,io + {Cf^^ + C'^)aLM,u> - i^lL{^)aLM,LO = Ilm^u,- (2.26) 

The collective mass -Bl is found to be 

r 3 

BL = m J drpeqJ2WLM,u\'^ = (2.27) 



The static stiffness coefficient Cj^^^ is derived from the first term on the right hand side of 
Eq. (p.l6|) and is given by |0] 

where bs and be are the surface energy and Coulomb energy coefficients appearing in the 
nuclear mass formula, respectively. This definition coincides with the one for the stiffness 
coefficient in the traditional liquid drop model for the nucleus. We point out, that the 
nucleon-nucleon interaction, manifested at the starting equations ( p.l| ) and (|2.2|), is presented 
in Eq. ( |2.26| ) only implicitly through the phenomenological stiffness coefficient Cj^^\ 

At finite frequencies, the distortion of the Fermi surface causes an additional contribution, 

in Eq. (|2.26| ), to the stiffness coefficient. It is found as 



C, ^ C'^iu) = Im (y^^) / rfrP,,Ai^''V,a2^,,, (2.29) 



where 



-r-(LM) ^ ^ 



Using a^M,!/ from Eq. ( |2.25| ), we obtain for the integral in Eq. ( [^.29D 



j drp^.Al'lf^V.alj,,^, = dL Po (2.30) 
where the information about the multipolarity is in 

Thus we find 

Ciiu) = dLlm (—^) Per (2.31) 

The proportionahty to {ujtY-, for small values of this product, explains, why such a correction 
does not appear in the hydrodynamic limit. 

For the friction coefficient 7l(co') in Eq. ( p.26|) we obtain 



7l(^) = j drAll^'^V^alj,,^, = rf^r^(^). (2.32) 

Both, C£ and ■jl depend implicitly on the temperature via the dependence of the collision 
time r and of Rq Peg on the T. 



The random force fLM,uj in Eq. (|2.26|) is related to the random pressure tensor Sj^^^^j by 



fLM,oj = -ni J dr SufM,oj'^ f^alj^j ,^. (2.33) 



Using Eqs. ( p.20|) and ( |2.23|) we obtain the spectral correlation function {fLM)t of the 
random force fiMit): 

'[h^ = 2Eiuj,T)viuj) I rfrAl^*'VXA/,. = 2 Eiuj,T) ^^H. (2.34) 

The basic property of the random variable y, in Eq. y = yjjl = transfers to both. 



the random pressure tensor, s^^^^ = 0, and the random force, /iAf.w = 0. 



Finally, using Eqs. ( p.26| ) and ( |2.34|) we can derive the spectral correlation function 



[01)1 for shape fluctuations 



2E{u;,T)-fL{uj) 



The result Eq. (2.35) is of similar form as one would obtain starting from hydrodynamics 
with the difference that here both the viscosity and the eigenfrequency ujl of the under- 
damped oscillator, 



+ (2.36) 



depend on the frequency themselves. This u dependence is due to the deformation of the 
Fermi surface and the non-Markovian collision term. 



III. NON-MARKOVIAN COLLISION TERM 

We have already pointed out, that memory effects in the collision term causes the collision 
time r, introduced in ( |2.14| ), to depend on frequency. In this section, we want to discuss 
this dependence in more detail. 

Below, we will take into account the main contributions to the relaxation time r, one 
arising from interparticle collisions (two-body dissipation with relaxation time T2) and the 
other from collisions of nucleons with the moving nuclear surface (one-body dissipation with 
relaxation time ri). Thus, 

111 , , 

- = - + -. 3.1 

The introduction of one-body dissipation reflects the peculiarities of our consideration: 
1) The presence of the sharp edge in the finite nuclear Fermi-liquid drop leads to a renormal- 
ization of the collisional integral I[6n] which is traditionally taken from Fermi-liquid theory 
of an infinite system; 2) The restriction to Fermi-surface distortions of multipolarity / < 2, 
see eq. (|2.2|), does not allow us to take into consideration Landau damping. Thus the as- 



sociated fragmentation width of collective excitations is missing. Introduction of one-body 



dissipation provides a phenomenological description of both contributions to the relaxation 
time. 

The Ti is related to the partial width T^^^ describing the damping of the collective state 
due to one-body dissipation, 

n = -J-,. (3.2) 



J- L 



For the modified one-body dissipation [^|[ the width r^^"* is given by 

V^l^ = ^T^ [L -\f L p^vpn\\ (3.3) 

where A is a free parameter which we will take from a fit to the experimental data of the 
width of the multipole giant resonance (MGR). For the case of finite temperatures, the Fermi 
velocity vp Eq. (p.3|) has to be replaced by the temperature dependent value, see Ref. 
0. The corresponding temperature dependent correction to the one-body relaxation time 
Ti, Eq. ( |3.2| ), is small at T << e^? and we will be neglected in this work. 

Let us now discuss the contribution of two particle collisions. One purpose of the present 
paper is to study modifications in the expressions for transport coefficients and fluctuations 
due to memory effects in the collision integral. These memory effects are realized in the 
dependence of not only on the distribution 6n at a given time but also on the value 
of during earlier times. The weight with which the distribution bn at previous times t' 
contributes to the value of at a given time t is given by the kernel of the convolution 
integral representing /[5n]. In the present case, this kernel depends on the difference t — t' 
only, since we have linearized the collision integral with respect to equilibrium. 

The inverse collision time 1/t2, Eq. (|2.14|) , is the Fourier transform of the kernel men- 



tioned and thus depends on the frequency. As a consequence of causality, l/r2 is actually 
a complex function of cj, whose real and imaginary part are related by the Kramers-Kronig 
relations. Thus, inclusion of memory effects in the collision integral /^^ in Eq. ( |2.5| ) implies 
a consistent change of both dissipative and conservative forces in the equations of motion 
(|11) and 
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A more detailed analysis of this question will be the subject of further investigation. 
Here, we only want to mention that we adopt phenomenological values for the parameters in 
our approach. Therefore, for the frequency range in which we are interested in the present 
case, it seems possible to neglect the frequency dependent corrections to the conservative 
forces in Eq. (^^.261). 



This approximation corresponds to keeping the real part of l/r2 only, the inverse of 
which will, in the reminder, be called collisional relaxation time and be labeled by the same 
symbol, T2- For frequencies small compared to the Fermi energy, one finds 

"^ = "^("'^)=(M^TCT^ (3-4) 



where (3 and C, are constant. For infinite nuclear matter (3 can be related 10| to the differ- 



ential cross section for the scattering of two particles. The differences between the different 



estimates of the parameter j3 |Tl|-|T5| are rather large: 

/3 = 2.4 ^ 19.8Me\/. 
For C we will adopt the value of C = 4 vr^ [|],|[ . 

IV. ZERO AND FIRST SOUND LIMIT 



Equation ( 2.35 ) is valid for arbitrary collision times r and thus describes both the zero 



and the first sound limit as well as the intermediate cases. From it one can obtain the 
leading order terms in the different limits mentioned. 

(1) First sound limit: ujt 0, T » hu 
The contribution from the dynamic distortion of the Fermi surface can be neglected in this 
case and we have from Eq. ( p. 31 ), 



^ 0. (4.1) 

The eigenfrequencies ul of the shape oscillations are determined here by the usual liquid 
drop model as 
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^ ^0^)13, . (4.2) 

In the high temperature regime, the Fermi hquid viscosity Tj^uj), Eq. approaches the 

classical expression 

V = ^Popk(O), (4.3) 
5 

where pp is the Fermi momentum and r(0) = t{uj = 0). The spectral correlation function 



(/l)J of the random force can be found from Eqs. ( |2.34| ), ( p.21| ) and (|2.32|) 



ih)l = 2^L{0)T. (4.4) 

This correlation function is independent of u, i.e., it corresponds to a white noise. 

(2) Zero sound regime: ujt —>■ oo, T << huj 
The contribution to the stiffness coefficient from the dynamic distortion of the Fermi surface 
is now given by (see Eq. ( p.31| )) 



Ciiiu)^C'L = dLPe,. (4.5) 

This expression coincides with the analogous one from [Q. We note that, in a cold Fermi 
liquid drop, provides the main contribution to the stiffness coefficient. The viscosity 
coefficient ri{uj), Eq. ( p.l8|) , can be approximated in this limit by 



r]{uj) = (PejKo) [l + ^ (C + ^)] , /to = 4 vr^ /? ^ c. (4.6) 



The random force spectral correlation function (/l)J is obtained from Eqs. (|2.34| ), ( p.32| ) 



and (ICTI ) to be 



{h)l = ^^7L, (4.7) 

where 

7L = dL Peg/ Ho (4.8) 

does not depend on u. The spectral correlation function ( ^.7| ) now corresponds to a blue 
noise. 
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We recall that the quantum-mechanical zero-sound regime Eq. (|4.6|) was obtained from 
the classical approach. It is due to the fact that the quantum fluctuations have been incor- 
porated into the correlation function (p.20| ) through the factor E{uj,T), Eq. ( p.21| ), see also 
Refs. @,|]. 

V. RESULTS AND DISCUSSION 

For the numerical calculations in this work we adopt the value of ro = 1.12 fm and assume 
a temperature dependence of the surface and Coulomb parameters in the liquid drop stiffness 
coefficient Cf of Eq. (g;^), namely [[ 



bs = 17.2 (^^f'ueV, be = 0.7(1 -xcT^)MeV, (5.1) 



-2 



where the parameter xc was chosen as xc = 0.76 ■ 10 ^ MeV and Tc = 18MeV is 



taken as the critical temperature for infinite nuclear Fermi-liquid ||T^. Using Eq. ( pTT 



one can find a critical temperature T^^^^ where the liquid drop contribution Cj^^^ to the 
stiffness coefficient vanishes: 

c^r^=c^r^(T)|^^^(..)=o. (5.2) 

in 

For the parameters used in the present work one obtains T^^n^^ = 7.72 MeV for quadrupole 
deformation, L = 2, in "^^^Pb. For temperatures T > T^^^^ there exists always an eigenfre- 
quency with a positive imaginary part giving rise to an exponentially growing deformation. 
This eigenfrequency, along with possible others, are solutions of the following secular 



equation, see Eq. (|2.26|) 



-u''Bl + (Cf + C'l{uj)) -iu^l{uj) = 0. (5.3) 

The transport coefficients C'j^{uj) and 7l(ci;) are cu-dependent because of the memory effects. 
To solve Eq. ( |5.3| ), both coefficients have to be defined in the complex cu-plane through 
analytical continuation of the corresponding expressions ( |2.31D and ( ^^.321 ). In Fig. 1 we 
show the real (Reu;) and imaginary (Imtu) parts of the eigenfrequencies, obtained by solving 
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Eq. (|5.3| ) as functions of the temperature for the nucleus with A = 208, using ( = in'^, Eq. 
( p.4| ), and three different values of l3. The value of A, see Eq. ( |3.3| ), was fixed by calculating 
the imaginary part of the eigenfrequency from Eq. ( p. 3D at T = and comparing the result 
with the experimental value of the half- width Tgqr/'2^ of the giant quadrupole resonance 
in a cold nucleus "^^^Ph [|T|]. We find X/tq = 3.011 for p = 3.4 MeV and X/vq = 3.551 for 
P = 9.2 MeV. The value of (3 = 1.26 MeV is the lower limit for (3 in the sense that, below 
this value, no positive value for A exists such that Imu equals the experimental value for 

rGQ/?/2. 

For each j3 there are three solutions to Eq. (5.3). One of them, ilmuin, is purely 
imaginary and lies, for T > T^^^\ in the upper half plane of the complex u (unstable 
mode). The two other solutions lie symmetric (for small enough T) with respect to the 
imaginary axis, iRecj — ihivuj, and meet at a temperature Tum- For example, we find 
Tiim = 9.5 MeV for f3 = 3.4 MeV. For larger temperatures, T > Tum, both solutions are 
imaginary and only the ones nearer to the real axis are shown in Fig. 1. Thus, for T > Tu^ 
the eigenexcitations exist as overdamped modes only. In Fig. 2 we have plotted Tu^ as 
function of the coUisional parameter f3. For (3 larger than about 2.3 MeV, the transition 
between underdamped and overdamped modes occurs in the presence of an unstable mode. 



In Fig. 3 we have plotted the spectral correlation function {aL)t as obtained from Eq. 
( p.35| ) for the two temperatures T = 1 MeV and T = 9 MeV. The different curves show the 



sensitivity of (cei)'^ to the parameter /5, Eq. ( |3.4|) . For low temperature we can observe a 
well defined maximum which corresponds to the GMR excitation (zero-sound regime). An 



increase of T leads to a shift of the maximum of (c^l)^ to lower frequencies and to an increase 
in the width. The shape of the curves near the zero-sound maximum is a non-Lorentzian 
one and depends, in particular, on the retardation effects in the friction coefficient, Eq. 
( p.32| ), and, consequently, on the parameters /3 and ( in the relaxation time, Eq. (^.41). 
Increasing the temperature we do not find a first sound peak centered at finite frequency 
for temperatures below T = t/^^^ and for realistic values of (3, see Fig. 3. The strong 
increase at low frequencies is due to the purely imaginary eigenfrequency. Thus, there forms 
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no resonance structure of the spectral correlation function (aL)S in this region of uo. 

Moreover, our numerical calculations, with realistic values of the nuclear parameters, do 
not show any transition from the zero-sound regime at low temperatures to the first-sound 
regime at high temperatures. For the finite Fermi-liquid drop, in the first-sound regime the 
real part of the eigenfrequency is mainly determined by the liquid drop stiffness coefficient 
Cj^^\ i.e., without the additional contribution C'^ from the Fermi-surface distortion effect. 
As it can be seen from Figs. 1 and 2, even if the limiting temperature is well below 
the temperature t'>^^\ the first sound regime does not appear because the eigenmotion is 
overdamped at < T < T^n^\ 

For a large enough value of /3, i.e. P > 2.3 MeV, there is, in principle, a possibility 
for a resonance structure of (q;l)^ at temperatures T > t/^^'' from the pure Fermi-surface 
vibration in the momentum space. For these values of (3 there exists a temperature region 
^/rf^"* < T < Tiim where Ci^^\t) < and Cl{ujl) > 0, simultaneously. This implies the 
existence, in this high temperature region, of a particular eigenmode of the Fermi liquid 
drop where the restoring force is exclusively due to the dynamical Fermi-surface distortion. 
Unfortunately, we find this eigenmode to be damped too strongly to cause a visible peak 
in the spectrum of the fluctuations for values of (3 < 25 MeV. As shown in Fig. 3, (q;l)^ 
develops only a weak shoulder for the value (3 = 9.2 MeV at T = 9 MeV. 



VI. SUMMARY AND CONCLUSIONS 

Starting from the collisional Landau- Vlasov kinetic equation with a random force, we 
have derived the Langevin equation for the shape fluctuations in the framework of the 
nuclear Fermi-liquid drop model. The main features of these fluctuations are due to the 
Fermi-surface distortion effects. We have obtained the random-force correlation function 
( p.34| ) for the general case where retardation and memory effects in the pressure tensor 
and in the relaxation time ( |3.4| ) are taken into account. It is important to account 
for these memory effects in order to obtain a correct description of the transition in Fermi- 
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liquid from the zero-sound regime at low temperature to the first-sound regime at high 
temperatures. We have found, however, that in a finite nuclear Fermi-liquid the first-sound 
regime is not reached. Instead, in the region of reasonable values for our parameters, the 
eigenmotion becomes overdamped at a temperature Tum. For Tum below the temperature 
t/^^-* the reappearance of underdamped motion is suppressed since the temperature for 
which C'^ = is smaller but close to Tu^. Note also that at temperature t/^^^ one of the 
modes becomes unstable. 

Our approach to the shape fluctuations is essentially classical. However, due to the 
quantum version of the fluctuation dissipation theorem , basic quantum effects are taken 
into account. Thus the correlation functions ( p.20| ) and ( p.34| ) contain contributions from 
both quantum and thermal fluctuations. 

The effects of the dynamical distortion of the Fermi surface on the nuclear collective 
motion lead to peculiarities of the random-force correlation function which are absent in a 
classical system. The spectral correlation function ( p.34| ) is independent of uj and corresponds 
to a white noise in the first-sound regime at lot 0, whereas in the opposite zero-sound 
regime at ujt — > oo it corresponds to a blue noise 



The behaviour of the spectral correlation function {dLit at different temperatures reflects 
the above mentioned peculiarities of the random-force correlation function. The dependence 



of the shape of the curves (aL)^ on the retardation effects (cj-dependence) in the friction 
coefficient ( p.32|) and, consequently, on the parameters C, and (3 may serve as means to learn 
about the possible range of the latter values. We find also the existence of an oscillatory mode 
due to the dynamical Fermi-surface distortion in the high temperature region T > tI^^\ 
where the usual liquid-drop stiffness coefficient C^"^^ disappears. The spectrum of the 
fluctuations, however, is affected by this mode only for collision frequencies which seem too 
high for stable nuclei. 
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APPENDIX 



In this appendix we are going to determine the correlation function of the projected 
random force yu^{j,t) Eq. (|27 



To calculate the correlation function of the random force y = yp{r,t), appearing in Eq. 

we follow the arguments of Abrikosov and Khalatnikov However, we generalize the 
treatment with respect to the following points: we allow for a collision term which is (i) non 
Markovian and (ii) nonlocal in space: 

= J ^^^^^J{T,p,t-r\p\t')6np,{T',t'), (A2) 

where 

(5np(r, t) = (5np(r, t) - -gf-^f^P^ (-4-3) 

Up is the equilibrium distribution function, 71° = 1 — n^, 6ep = J2p'r' /pp' (J" ~ r')Snpf (r', t) and 
/pp'(r — I"') is the quasiparticle interaction energy These generalizations do not change 
the expression for the rate of change of the entropy S as compared to [||]. In the present 
notation we write 

drdpdr'dp' 6np{r,t) /-^^^^^ „f\ 



^=-J (2nh)^ nOpffp ^(P-P)^(r-r)-^/ppKr-r) 

x[l[6n]p,y.t'+yp'ir',t')). (A.4) 



In Eq. (A.4) we corrected a typing error in the second term of the right hand side of Eq. 
(11.9) of Ref. 

To apply the general theory of fluctuations we need to write 5" in the form 
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In the present application the index i stands for the pair of vectors p, r and the sum has to 
be replaced by the integral over phase space. Defining 



Xi{t) = I[5n\r^^p^,t + ypA^i,t), {A.6) 

the expression (A. 3) for S takes the form (A. 5) if we set the generalized force Xi{t) = 
X{ri, Pi] t) as 

"'Pi "'Pi ^^tp^ "pi"pi 

In the next step, we have to express the first term on the right hand side of Eq. (A.6) in 
terms of the X(r, p; t), Eq. (A.7). With the help of Eq. (A.2) we find: 

I[Sn]r,p;t = Ir,p;t * {uy^X {v , p; t)) . {A.8) 

Comparing the general form of Xi{t), generalized to non Markovian processes 

= -Y.i^j{t,t')Xj{t')+y,{t) (A9) 

with Eqs. (A.6) and (A.8) we can easily determine the coefficients '~fij{t,t') 

7(r, p, t; r', p', t') = J{v, p, t; r', p', t')<rr° . (A. 10) 

If we assume that J depends on the difference t — t' only, we obtain for the correlation 
function of the random force 

duo 



;[l/p(r,t);yp,(r',t)]+ = / ^e-(*-*')i?(.;, T) 



X [J(r, p; r', p'; oj)jr^ + J{j\ p'; r, p; -uo)-^i (All) 



where E{uj,T) is determined by Eq. ( p. 211 ). This expression is valid for a quite general 



collision term. In the present paper we do not need the full correlation function, but only 
its projection (A.l). For the latter we find 
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x[J^(r,p;r',p';^)-^ + J(r',p';r,p; -uj)j^]pWt,'- 



der. 



{A.12) 



The collision integral used in the text is assumed to have the properties ( |2.14] ) with ( p.SI) . 
From the latter we find the kernel J in (A. 2) to fulfil the relation: 

[ dip' , , ^, , , .dnl 1 



whereas 



/ 



J/(r ,P ;r,p;u;)T^ = 



der. 



due to conservation of particle number in collisions. With the help of 



dp 



■—Wp' 



1 - 2^ 

dp P^, dn° 



one finds 



^[y•yt^ir,t)■,y^,^>{r',t')]_ 



|^e-M-')^(,,T)[J- + -l- 
27r t(uj) Ti—uj) 



dp r r r 2 



Here 



i2'Khflbrn?^ dep 

dp p^ 



P. 



eq- 



{2T:hf Ihm? ' dtp 
In Fourier space the correlation function is given from (A. 15) as 



(A 14) 



(A15) 



(A16) 



-[?/^^(r,a;);?/^v(r',^')] + 



Pe 2 
T\UO ) 6 

and the one for the random pressure tensor ( ^.19 ) reads 



(A 17) 



ATi5{uj + uj')5{y - v')E{uj,T) 



Peg. 

1 + [utV m'^ J 



With the definition of the viscosity coefficient (|2.18| ) we obtain formula ( p. 201) . 
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FIGURES 

FIG. 1. The dependence of the real, Rea;, and imaginary, lmu> parts of the eigenfrequency uj 
from Eq. (5.3) on the tcmpcratmc T for a nucleus A = 208 for the quadrupolc vibrations L = 2; the 
relaxation time parameters (5 = 1.26 MeV (dashed lines), 3.4 MeV (solid lines) and /? = 9.2 MeV 
(dotted lines) are shown as a label to the curves; the parameter ^ is chosen as C = 47r^. 

FIG. 2. The limiting temperature Tujn as a function of the relaxation parameter /?; using 
C = 47r^. The doted line is the critical temperature T^^^^ where the liquid drop stiffness coefficient 
ci^^'* vanishes. 

FIG. 3. The spectral correlation function {ai)t for two temperatures: T = IMeV and 
T = 9 MeV. The calculations were performed for ^ = 47r^ and the values (3 = 3.4 MeV and 
/? = 9.2MeV. 
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